LSND, SN1987A, and CPT ViolationEI 
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We point out that neutrino events observed at Kamiokande and 1MB from SN1987A disfavor the neutrino 
oscillation parameters preferred by the LSND experiment. For Am 2 > (the light side), the electron neutrinos 
from the neutronization burst would be lost, while the first event at Kamiokande is quite likely to be due to an 
electron neutrino. For Am 2 < (the dark side), the average energy of the dominantly u e events is already 
lower than the theoretical expectations, which would get aggravated by a complete conversion from to z? e . 
If taken seriously, the LSND data are disfavored independent of the existence of a sterile neutrino. A possible 
remedy is CPT violation, which allows different mass spectra for neutrinos and anti-neutrinos and hence can 
accommodate atmospheric, solar and LSND data without a sterile neutrino. If this is the case, Mini-BooNE 
must run in v rather than the planned v mode to test the LSND signal. We speculate on a possible origin of CPT 
violation. 



The neutrino masses are strictly zero in the standard model, 
while recent strong evidence for oscillations in atmospheric 
neutrino data suggests a small but finite mass for neutrinos [[]]]. 
There are also weaker but compelling hints for oscillation in 
solar neutrino data [||]. Both of them rely on the "disappear- 
ance" of the neutrinos compared to theoretical expectations. 
On the other hand, there is a dedicated neutrino oscillation 
experiment, LSND, which reported the appearance of v e in 
the flux from the stopped /i + decay [§,01. They have also 
reported a hint for appearance in — > v e mode but the signif- 
icance is low [g]. It was reported that its significance became 
even lower in the final analysis W], and hence we will ignore 
this hint throughout this letter. It is therefore an important 
question if all three indications for neutrino oscillation would 
fit together. 

The observation of neutrinos 

from SN1987A at Kamiokande and 1MB marked the birth of 
neutrino astronomy, and confirmed the standard core collapse 
model of Type-II supernovae. Detailed comparisons of data 
and theory put constraints on neutrino oscillation parameters 
(see for a review and reference therein). It is the aim of 
this letter to reexamine the constraints from SN1987A neu- 
trino data with a particular focus on the oscillation parameters 
preferred by the LSND experiment. 

There are basically three types of constraints one can draw 
from the SN1987A data. The first constraint comes from the 
energy spectrum of observed events, which are believed to be 
dominated by v e events. Because of different reaction rates 
in the proto-neutron star core, one expects a temperature hi- 
erarchy T Ve < T Pc < T V{Afi)A ^p T . Their average energies 
are expected to be 10-12 MeV, 14-17 MeV, and 24-27 MeV, 
respectively. The observed energy spectrum at Kamiokande 
indicates that the temperature of v e was somewhat on the 
low side of the theoretical expectations, with an average en- 
ergy of 7-14 MeV [^]. If there is an efficient conversion be- 
tween v e and or D T , it would increase the energies of the 



D e -induced events, aggravating the tension between data and 
theory. Therefore, oscillation parameters that would lead to 
such an efficient conversion are disfavored [QH. The MSW 
effect via the resonance occurs when P e is heavier than 
for small mixing angles as suggested by the LSND data, be- 
cause the matter effect due to the charged-current interaction 
with the electrons would bring the instantaneous eigenvalue of 
the Hamiltonian of S'e state lower and hence can cause level 
crossing. For Am 2 = 0.1-1 eV 2 , as suggested by the LSND 
data, the conversion is essentially complete and therefore the 
SN1987A data disfavor such parameters. We use the require- 
ment P osc < 0.35 by Smirnov, Spergel and Bahcall based on 
this argument [Q]. However, the constraint had not been stud- 
ied on the dark side tan 2 6 > 1 of the parameter space [||] to 
the best of our knowledge Jl0| ] . The density profile was taken 
from Jll|] with an empirical approximation 



N e {r) 



IQ 10 N A 

O.lA^(r/r )- 3 



're 



(1) 



where Na is the Avogadro number, and we took E v ~ 
25 MeV. The resulting constraint on the oscillation parameter 
space is shown in Fig. [I] on the dark side. The wiggles around 
Am 2 ~ 10~ 5 eV 2 are due to the Earth matter effect. Because 
the Large Magellanic Cloud is seen on the southern sky while 
both Kamiokande and 1MB detectors reside on the northern 
hemisphere, the neutrinos from SN1987A had passed through 
the Earth, causing regeneration of v e and hence making the 
constraint weaker. We approximate the effect using a constant 
density N e ~ 3N A and R ~ 10 4 km. 
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FIG. 1. Constraints on the two-flavor oscillation parameter space 
from SN1987A neutrino data. The shaded MSW triangle on the light 
side (tan 8 < 1) is disfavored by the neutronization v e burst, while 
that on the dark side (tan 2 > 1) by the energy spectrum of the 
v £ induced events. The constraint in order to preserve the nuclear 
r-process excludes the region above the curve [|l4[. LSND preferred 
region is shown at 90 and 99% CL. KARMEN2, Bugey constraints 
at 90% CL are taken from [jlfj and exclude the regions above the 
curves. The currently preferred regions (95% and 99% CL) from the 
solar neutrino data W% are shown for comparison. 

The second constraint comes from the very first (and possi- 
bly the second) event at Kamiokande. At the time of core col- 
lapse, most of the protons in the iron core of the progenitor are 
converted to neutrons to overcome the Coulomb repulsion, re- 
leasing electron neutrinos. This is called the neutronization or 
deleptonization burst. Near thermal radiation of all species of 
neutrinos used in the first constraint appear only about a hun- 
dred milliseconds after the neutronization burst. The electron 
neutrinos dominantly scatter elastically with electrons in wa- 
ter, and produce highly forward peaked electrons, while T> e ab- 
sorption on proton produces a nearly isotropic ally distributed 
positrons. Indeed the very first event at Kamiokande points 
beautifully back at the SN1987A and is completely consistent 
with this interpretation. The expected event rate of v e is, how- 
ever, about 0.025 at Kamiokande [|l2[| and hence the observa- 
tion is thanks to an upward statistical fluctuation. If there is 
an efficient conversion between v e and or v T , the expected 
neutral-current event rates due to would be about 1/7 of 
that due to the v e events which have both neutral-current and 
charged-current amplitudes, and hence the observation of one 
event would be highly unlikely. Therefore, oscillation param- 
eters that would lead to such an efficient conversion are dis- 
favored [13 nj]. We require P osc < 0.90, so that the obser- 
vation of one event is possible within 99% CL. The MSW 
effect via the resonance occurs when v e is lighter than 
for small mixing angles as suggested by the LSND data, be- 
cause the matter effect due to the charged-current interaction 
with the electrons would bring the instantaneous eigenvalue of 
the Hamiltonian of v & state higher and hence can cause level 
crossing. For Am 2 = 0.1-1 eV 2 suggested by the LSND 



data, the conversion is essentially complete and therefore the 
SN1987A data disfavor such parameters. The constraint on 
the oscillation parameter space is shown in Fig. |l] on the left 
half tan 2 6 < 1 (the light side). 

The third constraint is based on the assumption that the ex- 
panding envelope driven by thermal neutrino wind of explod- 
ing supernova is the site of nuclear r-process, synthesizing 
elements beyond iron. If the v e is lighter than v^, there may 
be an efficient conversion between the two, and the v e wind 
would have the temperature of v^, i.e. higher than what it nor- 
mally is. The higher temperature of v e would have a higher 
cross section to convert neutrons to protons, where protons 
would end up mostly in 4 He and would not participate in pro- 
ducing neutron-rich nuclei required in the nuclear r-process. 
This consideration places a constraint at higher values of Am 2 
Jl4|Jl5| ]. The constraint derived in jl4} | is shown in Fig. [jj 

For a comparison, we also show the preferred regions of the 
parameter space from the solar neutrino data in Fig. [j], taken 
from [jlj]. 

The important point is that the oscillation parameters pre- 
ferred by the LSND data, which could be on either sides of the 
parameter space, are both disfavored by the SN1987A neu- 
trino data, even though the difficulties in theory of supernova 
explosion and low statistics in the data do not allow us to draw 
a definite conclusion. 

So far, our analysis has been within the two-flavor mixing 
scheme. However we know there are three light active neu- 
trinos, and it has been argued that we may need even a sterile 
neutrino state to explain LSND, atmoshperic, and solar neu- 
trino data by neutrino oscillations. Note that our result does 
not depend on other oscillation effects, in particular whether 
there exists a sterile neutrino or not. Let us consider the case 
where all current indications for neutrino oscillations, LSND, 
atmospheric, and solar, are correct and hence there is one ster- 
ile state. Because there are three independent Am 2 , there are 
3! = 6 ways to order them. Four of them are so-called 3+1 
models, where one state is separated by Atol SND while other 
three are close to each other separated only by Am 
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These models used to be disfavored by the com- 



bination of CDHS, CCFR, Bugey, and atmospheric neutrino 
data at SuperKamiokande [[18|], but recent reanalysis of the 
LSND data brought the preferred Am 2 and sin 2 29 smaller 
and there opened a small acceptable region in the parame- 
ter space jigp. In these models, the state v± widely sepa- 
rated from the rest is nearly pure u a , with small mixing of 
v e and v^. The LSND oscillation is explained by the prod- 
uct sin 2 20lsnd = 4|£/ e 4[/* 4 | 2 . Two other models are so- 
called 2+2 models, where two doublets, each responsible for 
atmospheric and solar neutrino oscillations, are separated by 
Am 2 SND . v e (Vfj) state is almost exclusively in the solar (at- 
mospheric) doublet. Recent SuperKamiokande data disfavors 
pure v s oscillation in both solar and atmospheric data, and 
therefore we need to put v s in both doublets [pC|]. Now we 
follow how the states evolve as the neutrinos exit the proto- 
neutron star core. In 3+1 models, either v e or v e crosses the 
v s state first, and the transition between these states is well ap- 
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proximated by the two-flavor mixing as studied above. There- 
fore either of them is nearly completely lost into v s . This 
strengthens the constraint from the neutronization burst, while 
the constraint from the v e spectrum is unchanged by the v s 
component (even though the overall normalization gets fur- 
ther suppressed). In 2+2 models, either v e or v e crosses the 
atmoshperic doublet first and are nearly completely converted. 
Therefore the constraints discussed in the two-flavor case are 
unaffected. 

A fair question to ask is how robust these constraints are. 
As for the temperature difference used in the first constraint, 
the issue had been raised if an additional process, such as 
vNN — > vNN, may reduce the temperature differences be- 



tween v e and v^, v^, v T , v T [|21j, but no concrete estimates 
of the temperature had been given. This issue can be settled 
only by more detailed numerical simulations and/or a future 
observation of supernova neutrino bursts. For instance, Su- 
perKamiokande, SNO, Borexino, and KamLAND can detect 
v e via the charged-current reaction, while v e can also be de- 
tected at SNO via the charged-current, and all neutrino species 
at SNO, Borexino, and KamLAND via the neutral-current re- 
action (see [ j22| ] for a recent review on the experimental as- 
pects). Then we can test if there is a significant temperature 
difference among different event categories. The interpreta- 
tion that the first Kamiokande event, produced at the angle 
18° ± 18° in the forward direction, is due to the elastic v e e 
scattering of v e from the neutronization burst is also subject 
to a criticism. The expected event rate is low, and we are 
relying on a single event to place the constraint. The proba- 
bility that this event is due to an isotropically distributed v e 
event is about 3% [jO]. We find the fact that this event was 
the first quite suggestive of being a v e event. Again a future 
detection of supernova neutrinos would settle this issue. Be- 
cause of these possible criticisms, we cannot make a definite 
claim that the LSND data is incompatible with the SN1987A 
neutrino events. We can only say that the LSND preferred re- 
gion is disfavored by SN1987A data based on the assumptions 
made above. 

For the rest of the letter, we take the above constraints 
seriously, and we discuss how we may accommodate the 
LSND data despite the constraints. The only way to evade 
the SN1987A constraints is to assume v e is heavier than iv T 
while v e lighter than v^. T . Such a mass spectrum obviously 
violates CPT, but we do not see any other alternatives as long 
as we take the SN1987A constraints seriously. Once CPT is 
violated, in principle one may also consider the violation of 
Lorentz invariance, which will be discussed elsewhere. For a 
phenomenological exercise, we consider different mass spec- 
tra for neutrinos and anti-neutrinos and keep Lorentz invari- 
ance. The question is if we can accommodate atmospheric, 
solar, and LSND data within the SN1987A constraints. 

The key to this question is that the solar neutrino data probe 
only v e , but not v e , while the LSND data only v e , but not v e 
J23]]. On the neutrino spectrum, unless c7 e3 element is ex- 
tremely small, v e cannot be below the atmospheric neutrino 
mass gap because it would cause a loss in the neutronization 



burst. Therefore, the solar mass gap must be above the atmo- 
spheric mass gap. On the anti-neutrino spectrum, we need v e 
the lightest, and v^, v T about 1 eV 2 above v e . The splitting be- 
tween two mass eigenstates which are dominantly v^, v T must 
be relevant to atmospheric neutrino oscillation. The mass 
spectra are depicted in Fig. ^. It is interesting to note that the 
combination of the cosmic microwave background anisotropy 
from Planck satellite and Lyman-a power spectrum will be 
able to exclude the neutrino mass down to 0.29 eV at 2a level 

Different mass spectra between neutrinos and anti- 
neutrinos will affect future neutrino oscillation experiments. 
The most important consequence is for the Mini-BooNE ex- 
periment, which is supposed to put a final word on the LSND 
signal. They will run primarily in the mode, unfortunately, 
which would not exhibit the LSND oscillation. They do have 
a capability to run in the mode, however, and this mode 
must be used to test the LSND evidence. 
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FIG. 2. Possible mass spectra of neutrinos and anti-neutrinos con- 
sistent with solar, atmospheric, LSND data and the SN1987A con- 
straints. 



Brief comments on the possible origin of CPT violation are 
in order. First of all, it has been often argued that the small 
mass of neutrinos could well be originated from Planck-scale 
physics. Even though Yukawa couplings suppressed by the 
nominal Planck h ~ v/Mpi would be too small compared 
to the required mass spectra above, the "Planck scale" can 
well be much lower, even down to the TeV scale as has been 
discussed intensively lately $3j. Therefore, it is quite pos- 
sible that the small neutrino masses probe quantum gravity 
physics. It is then also conceivable that the possible violation 
of CPT from quantum gravitational physics appears most ev- 
idently in the neutrino mass spectra but not elsewhere. For 
instance, non-commutative geometry violates Lorentz invari- 
ance at short distances, producing possible seeds for CPT vi- 
olation [ p6[ l ( see > however, [p7[]). It is easy to write down 
Hamiltonian with different massees for neutrinos and anti- 
neutrinos in momentum space, but it is non-local in the co- 
ordinate space. See Refs. [ p8[ ] for recent discussions on 
other stringy or quantum-gravitational origin of CPT viola- 
tion. Even though this discussion is highly speculative, we 
hope that our work provokes more intensive discussions on 
the possible origin of CPT violation. 
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In summary, we discussed the SN1987A constraints on the 
neutrino oscillation parameters, and found that the parameters 
preferred by the LSND data are disfavored by the SN1987A 
data on both sides of the parameter space. If we take these 
constraints seriously, the only way to make the LSND data 
compatible is to allow different mass spectra for neutrinos and 
anti-neutrinos, and hence CPT violation. 
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